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Tetraquarks are studied with an approximation which reduces the in- 
ternal degrees of freedom of the system, using both finite differences and 
scattering theory. The existence of bound states and resonances is inspected 
and the masses and widths of the found resonances are calculated. 
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1. Introduction 

For the qiq2q3Q4 system, there is evidence in Lattice QCD [Tl[2], at least 
for static quarks, that the hamiltonian, is 



H^J2^i + Vi234 (1) 



i=l 

with the potential being given by V1234 = ram{V^^*'''',V^^^^,V^^^*), 
where y^i^2 Vm^Ma are the two possible meson-meson potentials and 

yTetra = C^ + „^^^.^J_ + ^L^^„(ri, r2, Tg, (2) 

is the tetraquark sector potential. We call V1234 the triple flip-flop potential. 

In the tetraquark configuration, a single connected string is formed. 
This string is formed from five segments of string as can be seen in Fig. 
[1] (see [31 H] for a detailed explanation ). We have studied the colour field 
configurations for the tetraquark. The results are shown on Fig. [21 

It is hard to solve this hamiltonian, so we will instead solve a simplified 
potential model which contains the essential ingredients of the triple flip-flop 
potential, but depends on a smaller number of variables. 
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Fig. 1. Left: Triple flip-flop Potential potential. Right: Simplified potential, which 
depends only on p and r. 
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Fig. 2. Plot of Lagrangian density of a static tetraquark. ri is the intra-diquark 
distance and r2 the inter-diquark distance. 



We will do the following approximations. First we will consider that 
all the quarks and antiquarks have an equal mass, while neglecting the 
exchange effects. Secondly, we will consider only the confining (linear) part 
of the potential and that the quarks are non-relativistic. Finally, we will 
also reduce the number of variables of the model. 

2. A toy model for the potential 

We now simplify the triple flip-flop potential, in order to reduce it to 
a two-dimension problem, similar to the one of a cherry in a broken glass. 
To simplify the triple flip-flop potential, we assume that the quark-quark 
distance = 1*1—^2 and the antiquark-antiquark distance P34 = ra — r4 
are equal — P34 = P- Thus, we only need to consider the variable p and 

Y _ rrir2 _ ^^3+r4 model. 

In this way, only one meson-meson branch is relevant, and the sim- 
plified flip-flop potential is given by Vi?i?(r, p) = min(VMM5^)5 where 
^mm(^5P) = (j{2r) and VT{r^p) — a{r + VSp) This potential is depicted 
in the right side of Fig. [TJ 
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Fig. 3. Left: Semi-localized state, or resonance for Ir = 1. Center: Bound state for 
Ir = 3. Right: "Phase shifts" obtained from the finite differences. 



In this way, we obtain the Schrodinger equation, for this model: 



i^{r, p) = E^{r, p) ■ 



(3) 



and what we call the outgoing spherical wave method (useful to study res- 
onances) . 



3. Finite difference method 



Since the Hamiltonian, is invariant for independent rotations of p and 
r, we have two conserved angular momenta = r x and Lp = p x pp. 

So, the wavefunction can be written as = ^^^^^f^Yirmri^ri ^r)Yipmpi^p, ^p) 
and the Schrodinger equation® becomes 



2m V dr^ dp'^ 



u + 



U + 



-u + VpFix, p)u = Eu . (4) 



We now place the system in a lattice and discretize the second derivatives 
the usual way. Dirichlet boundary conditions are used at the origin, for the 
wavefunction to be regular and at the box boundary for convenience. 

We then evaluate the eigenstates of the system. Only a few of them are 
bound or semi-localized states. They can be identified by calculating v^(p^ 
and comparing it to the box size. 

In the following results, we consider only Ip — Q. For = 0, only 
continuum states are found. However, for 1^ > 0, we are able to found some 
semi-localized states, such as the one depicted on the left side of Fig. [3l On 
the center, we can see a true bound state found for 1^ = 3. 

We try to calculate the phase shifts using the finite differences results, 
by calculating the functions ipi{p) = / (i^r(/)*(r) ^^(r, p) where the are 

the eigenfunctions of the hamiltonian Hmm — "^^r + ^ar and, then by 
fitting the long range limit Ai sm{kip + 6i). We obtain the results in 

the right side of Fig. [3l As we can see there, the results become irregular 
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when a second channel is open, due the interference of the different channels, 
and so we can't use this method to calculate the phase shifts. 
For that, we had to resort to a scattering theory approach. 

4. Outgoing spherical wave method 

By expanding the wavefunction as ^^(p, r) = ^i^ii^^ p) we obtain the 
coupled channel Schrodinger equation 

- 1^ V^*, + V.j^j = {E- 6,)*, , (5) 

where Vij = / d^r^^ * {'^){Vff — 2crr)(/)j(r), are the eigenvalues of Hmm 
and jii — m. 

Considering the scattering from the channel i to the channel j, we have 

• 7 ikjr 

asymptotically "ifj 5ije^'^^^ ^ fij{r)^—^ ^ if the channel j is open, otherwise 
vanishing. The fij can be computed by considering the outgoing solutions 

of the Eq. © , = e'^^'5,j + Y.nm ^Ylm{0, We obtain [5] 

- + ^j^^^ = - ^j>j - Vjnji{kir)r . (6) 

Also we can obtain the phase shifts by studying the asymptotic behaviour 
of the solutions. 

We can describe the scattering process with four quantum numbers: The 
scattering angular momentum Zp, the confined angular momentum 1^ and 
the initial and final states radial number in the confined coordinate r, rii 
and rij. 

We now compute the phase shifts, in order to search for resonances in 
our simplified flip-flop model. 

The energies in the plots will be given in terms of £"0 = (/i^cr^/m)-'^/^, 
which is the energy scale of the problem. 

In Fig. Hlwe see the phase shifts for different values of from the radial 
ground state = 0. As can be seen, no resonance is formed for 1^ = 0, but 
for = 1 and Ij. — 2 one resonance is formed since the phase shift has a 
jump of TT. This is somewhat expected, since the centrifugal barrier has the 
effect of maintaining the diquarks separated, while increasing the stability 
of the system. We also see a bound state for 1^ — 3, which agrees with the 
results obtained with the finite differences method and presented in Fig. [3l 
The phase shifts for scattering for rir — 1 and different and for = and 
different values of are also presented. 

If we include the relativistic corrections for the reduced mass of the 
system, the fii becomes jii — rn+j^. Greater reduced masses could increase 
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(E-4mc')/E„ 



Fig. 4. Comparison of the phase shifts for = 0, 1, 2 and 3, with = (left) and 
rir = I (center). In the right, phase shifts, for = and different n^. 



1^ = 




Fig. 5. Phase shifts, including the reduced mass correction for = 0, 1 and 2. 



the binding of the system. Which is in fact what we observe in Fig. [5l For 
= we can see a resonance that is destroyed for masses as low as 2v^/c^, 
while for = 1 the resonance is always present. When 1^ — 2 for sufficiently 
low masses we can observe a bound state and even a resonance in addition 
to it. 

Finally we can compute the decay widths of the resonances by using the 
expression ^ — 

(H) wh^^ ^ is 7r/2. The results are given in Table [H 
The width is larger for large angular momentum and small quark mass. 

Using this results to estimate the physical scale of our decay widths, we 
obtain, in the Ir — 1 state, a decay width F/ 2^ 15MeV for light quarks 
(m 2^ 0.4 GeV). For the charm quark we obtain Fc ^ lOMeV and for the 
bottom quark, the estimated width is F5 2^ 6MeV. 

In the Ir — 2 state the estimated widths are F/ = 52MeV, Fc = 34MeV 
and F5 = 23MeV, while in the = 3 state, the values are F/ = 140MeV, 
Fc = 90MeV and F5 = 60MeV. 

5. Conclusion and outlook 

We have constructed a simplified model for the tetraquark. With this 
model, we have observed the formation of a resonance for non null and 
even a bound state for = 3. So, we conclude that a higher angular mo- 
mentum favours the formation of resonances and bound states, in agreement 
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{E -4mc^)/EQ 


r / Eo 


1 


6.116 


0.037 


2 


6.855 


0.131 


3 


7.462 


0.352 







{E - Amc^)lEQ 


r/ ^0 





1.0 


5.001 


0.039 




1.5 


5.096 


0.022 


1 


1.0 


5.659 


0.137 




3.0 


5.990 


0.075 




5.0 


6.053 


0.053 


2 


1.0 


6.194 


0.586 




3.0 


6.634 


0.209 




7.0 


6.777 


0.162 



Table 1. Left: Decay widths as a function of Right: Decay widths with the 
reduced mass correction. 



with Karliner and Lipkin Extrapolating this results to the real world, it 
seems plausible the existence of tetraquarks with high angular momentum. 

The formation of resonances and bound states is also favoured by low 
quark masses in our model. This happens since a low quark mass produces 
a higher binding energy which increasing the meson masses and so favouring 
the localization of resonances and bound states. 

We think the methods developed here are useful and could be futurely 
used in the complete tetraquark problem with a triple flip-flop potential. 

We thank George Rupp for useful discussions. We acknowledge the 
financial support of the FCT grants CFTP, CERN/FP/109327/2009 and 
CERN/FP/109307/2009. 
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